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ABSTRACT 

We perform collisionless A-body simulations to investigate whether binary mergers between 
rotationally-supported dwarfs can lead to the formation of dwarf spheroidal galaxies (dSphs). Our 
simulation campaign is based on a hybrid approach combining cosmological simulations and controlled 
numerical experiments. We select merger events from a Constrained Local UniversE (CLUES) simu- 
lation of the Local Group (LG) and record the properties of the interacting dwarf-sized halos. This 
information is subsequently used to seed controlled experiments of binary encounters between dwarf 
galaxies consisting of exponential stellar disks embedded in cosmologically-motivated dark matter ha- 
los. These simulations are designed to reproduce eight cosmological merger events, with initial masses 
of the interacting systems in the range ~ (5 — 60) x 10' Mq, occurring quite early in the history of 
the LG, more than 10 Gyr ago. We compute the properties of the merger remnants as a distant 
observer would and demonstrate that at least three of the simulated encounters produce systems with 
kinematic and structural properties akin to those of the classic dSphs in the LG. Tracing the history 
of the remnants in the cosmological simulation to z = 0, we find that two dSph-like objects remain 
isolated at distances > 800 kpc from either the Milky Way or M31. These systems constitute plausible 
counterparts of the remote dSphs Cetus and Tucana which reside in the LG outskirts, far from the 
tidal influence of the primary galaxies. We conclude that merging of rotationally-supported dwarfs 
represents a viable mechanism for the formation of dSphs in the LG and similar environments. 
Subject headings: galaxies: dwarf - galaxies: Local Group - galaxies: fundamental parameters - 
galaxies: kinematics and dynamics - methods: numerical 



1. INTRODUCTION 

The origin of dwarf spheroidal galaxies (dSphs) in the 
Local Group (LG) is a long-standing subject of debate. 
Thes e intriguing sy stems are the faintest galaxies known 
(e.g., iMateol fl998[ ) and they are be lieved to be highly 
dark matter (DM) domi nated (e.g.. lMateo1ll998t iLokasI 
200ft iWalkeret al.ll2009h . Moreove r dSphs are gas poor 
or co mpletely devoid of gas (e.g., iGrcevich fc Putmanl 
2009) and they are characterized by pressure-sup ported, 
spheroidal stellar components (e.g.. lMatedll998l) and a 
wide diversity of star formation histories (e.g., IGrebell 
2000) . While our understanding of dSphs has grown im- 
pressively in the past decade, a definitive model for their 
formation still remains elusive. 

One class of models attempts to explain the ori- 
gin of dSphs via various environmental mechanisms, 

including tidal an d ram pressu re stri pping (e.g. 

Einasto et all 1 1974: iFaber fc Lid 119831: 
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20011: iKravtsov et al.1 12004: iMaver et al.l 
Klimen towski et al.l | 200ft iK azantzid is et al.1 120 lTT ) and 
resonant stripping ()D'Onghia et al.l I2009D . In th is con- 
text, the tidal stirring model (jMaver et aT1l200lD posits 
the formation of dSphs via the tidal interactions between 
rotationally-supported dwarfs and Milky Way (MW)- 
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sized host galaxies. While this model naturally explains 
the tendency of dSphs to be concentrated near the dom- 
inant spirals, it is challenged by the presence of the iso- 
lated dSphs Cetus and Tucana in the LG outskirts. Such 
examples illustrate the need for complementary models 
concerning the origin of dSphs. 

In the current galaxy formation paradigm (e.g., 
I White fc Reeslll978T ). the quiescent cooling of gas within 
a virialized DM halo results in the formation of a 
rotationally-supported disk of stars. Thus, any scenario 
for dSph formation must incorporate physical processes 
for transforming initially rotationally-supported stellar 
systems to ones dominated by random motions. In- 
teractions and mergers between galaxies may constitute 
such a mechanism. Indeed, on larger scales, the tidal 
heating and violent relaxation associated with mergers 
of massive, disk galaxies effectively destroy the stellar 
disks creating a kinematically hot, pressure-sup ported 
spheroid that resembles an elliptical galaxy fe.g.. lBarnesl 
Il992h . It has also been demonstrated that encounters 
between d warf halos can lead to their very strong evolu - 
tion (e.g.JKnebe et al.l 120061: iKlimentowski et al.ll2010h . 
iKlimentowski et al.l (|2010f ) reported in their constrained 
DM cosmological simulation of the LG that a few per- 
cent of all surviving subhalos had undergone a substan- 
tial encounter with another dwarf halo in the past. Most 
of these events occurred early in the LG history and be- 
fore the dwarf halos were accreted and became satellites 
of either the MW or M31. Therefore, although the great 
majority of LG dwarf galaxies should not have partici- 
pated in mergers with other dwarfs in the past, at least 
some of them may have experienced such interactions. 

Motivated by these findings, we selected eight merger 
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events from the same LG simulation and re-simulated 
them at higher resolution, embedding stellar disks inside 
the dwarf DM halos. Our goal is to determine whether 
mergers of initially rotationally-supported dwarfs can 
produce systems with kinematic and structural proper- 
ties akin to those of the classic LG dSphs. 

2. METHODS 

The basis for our experiments is a constrained DM cos- 
mological simulation of the Local Universe which repro- 
duces th e main LG properties at z = 0. This simula- 
tion (see[Klimcntowski et al. 2010]) is part of the CLUES 
project 5 (jGottlober et al.ll2010[ ). We identify dwarf-sized 
halos participating in mergers and record their proper- 
ties by utilizing the halo catalogues of the cosmological 
simulation th at were generated with the AMIGA halo 
finder (AHF) (jKnollmann fc Knebdl2009l) . 

In practice, it is not feasible to simulate all possible 
encounters identified in the halo catalogues at sufficient 
numerical resolution. For our purposes, it would be suf- 
ficient to show that even a single representative merger 
between dwarf galaxies can lead to the formation of a 
dSph. To this end, we selected eight binary interactions, 
which we denote M1-M8, for subsequent re-simulation. 
These encounters corresponded to mass ratios of the in- 
teracting systems ranging approximately from 1:1 to 1:4. 
While these merger events were not chosen in any special 
way, they were associated with a significant mass transfer 
between the systems (> 80%) indicating the occurence 
of an actual merger, as opposed to a weak encounter or 
a fly-by. 

Table [T] contains the properties of the two interact- 
ing systems in each case. Column 2 lists the redshifts 
at which the mergers were identified in the cosmological 
simulation. All recorded encounters took place at very 
early cosmic times (z > 2.5), in accordance with the ex- 
pecta tion for the epoch of r apid growth of low-mass halos 
(e.g.. lWechsler et al J 120021) . We note that these mergers 
occurred before the dwarf halos have been accreted onto 
more massive structures. Indeed, subsequent to accre- 
tion, the merging probability of halos di minishes due to 
the l arge host velocity dispersion (e.g., IDe Lucia et al.l 
12001 . 

Columns 3-6 list the halo virial parameters (M v - lr , R V1I , 
Vvir) and concentrations c = R V i r /r s (where r s denotes 
the "scale radius" where the product p(r) r 2 reaches its 
maximum value). Colu mns 7-10 list the dimensionless 
spin parameters A (e.g., [Peebles 1969) and the direction 
of the angular momentum vector of the dwarf halos nor- 
malized to unity. 

We constructed numer ical realizations of dwarfs using 
the method described in lKa zantzidi s~et al.l (|2005). The 
models consisted of expon ential stellar di s ks em bedded in 
cosmologically-motivated iNavarro et al.l (|1996| ) DM ha- 
los whose properties match those of the cosmological ha- 
los in Table [TJ We parametrized the mass of the stellar 
disk as a fraction, fd, of the halo M v ; r and the disk thick- 
ness as Zd/Rd {zd and Rd denote the disk vertical scale 
height and radial scale length, respectively) . For simplic- 
ity, we fix the values of Zd/Rd and fd in all models. In 
our modeling, Rd is determ ined by M v j r , A, c, and fd via 
the semi- analytic model of I Mo et al.l (|1998|) for the for- 

5 http://www.clues-project.org 



mation of disks 6 . The resulting Rd are listed in column 
11 of Table [U 

We take the expectation that dwar f galaxies should 
be b orn as thick, puffy systems (e.g., IKaufmann et all 
120071) into account by conservatively adopting Zd/Rd = 
0.2. Moreover, we choose fd = 0.005. Such small 
value is consistent with that inferred for some of the 
faintest LG dwarf i rregular galax ies (dlrrs) such as L eo A 
(IBrown et al.ll2007l) and SagDIG dYoung fc Lolll997l) . As- 
suming an isothermal DM halo, the velocity dispersions 
of these systems (a ~ 10 kms -1 ) correspond to V V1I ~ 
14 kms^ 1 , comparable to those of our halos. Tiny stel- 
lar components are expected in such low-mass systems 
because both gas accretion/retention and star formation 
would be strongly suppressed by re-ionizat io n and su- 
perno vae driven outflows (e.g., Maver1 l2"010D . lOh et all 
(|201lD computed fd for a fairly large sample of dwarfs in 
the THINGS survey, obtaining values ~ 0.01, albeit for 
galaxies that would have been less affected by the afore- 
mentioned effects having rotational velocities that are a 
factor of 2 — 3 larger than those considered here. 

We rotated each dwarf model so that the direction of 
the angular momentum vector of its DM halo was identi- 
cal to that of the corresponding cosmological halo. In our 
modeling, the angular momentum of the disk is aligned 
with that of the host DM halo. This results in rela- 
tive orientations of the disk angular momenta that span 
a wide range (~ 10° — 120°). Lastly, we realized the 
merger simulations using the initial positions and veloci- 
ties of the corresponding cosmological halos with respect 
to their center of mass (Tabled)). 

All merger simulations were perfor med with the N- 
body code PKDGRAV (IStadell I200ll) . Each A-body 
dwarf model contained Nh = 3 x 10 5 DM particles and 
Nd = 5 x 10 4 disk particles. The gravitational softening 
was set to eu — 30 pc and td = 15 pc, respectively. Force 
resolution was adequate to resolve all scales of interest. 

3. RESULTS 

All merger products were allowed to reach equilibrium 
before any analysis was performed. We first determined 
the principal axes of the stellar components using the mo- 
ments of the inertia tensor. Subsequently, we produced 
2D maps of the stellar surface distribution and kinemat- 
ics along the three principal axes and "observed" the 
merger remnants as a distant observer would. 

Given that o bserved dSph s are supported by random 
motions (e.g., Mateo 1998), only remnants satisfying 
Kot/cr* < 1 (where V to t and cr» denote the stellar ro- 
tational velocity and the one-dimensional, central veloc- 
ity dispersion, respectively) may be regarded as dSphs. 
Using the kinematic maps of the stellar mean radial ve- 
locity and velocity dispersion, we computed Vrot/c* for 
all remnants. For each line-of-sight, the value of V TO t cor- 
responded to the maximum velocity found anywhere on 
the map. For the velocity dispersion cr» we adopted the 
central value. 

6 We note that the IMo et aLl 1119981 ) formalism may be inap- 
propriate for dwarf galaxies at high-redshift. However, given the 
lack of knowledge f or the structure of dwarf galaxies at early cosmic 
times, utilizing the Mo et al. ( 1998) model constitutes a reasonable 
alternative to assigning arbitrary values to the disk scale lengths 
of our dwarfs. 
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TABLE 1 

Initial Properties of Merging Dwarfs 





M v i r 


-Rvir 


Vvir 








Rd 


Merger z 


(1O 7 M ) 


(kpc) 


(km/s) c A 


Lx 


Ly 




(pc) 


(1) (2) 


(3) 


(4) 


(5) (6) (7) 


(8) 


(9) 


(10) 


(11) 



Ml 



3.02 



6.51 
5.99 



3.15 
3.03 



9.43 
9.22 



4.48 0.046 
3.08 0.032 



-0.34 
-0.21 



0.73 
-0.22 



-0.59 
-0.95 



119 
97 



M2 



3.82 



7.31 
7.17 



2.69 
2.67 



10.82 
10.75 



4.11 0.039 
1.92 0.033 



-0.18 
-0.70 



-0.41 
0.68 



-0.89 
0.21 



92 
103 



M3 



6.67 



5.78 
7.34 



1.56 
1.69 



12.61 
13.66 



1.97 
1.46 



0.054 
0.047 



-0.05 
-0.11 



-0.86 
-0.76 



-0.50 
-0.64 



90 
91 



Ml 



3.15 



58.37 
19.95 



6.24 
4.36 



20.05 
14.03 



1.34 
1.13 



0.032 
0.067 



0.74 
0.25 



0.64 
0.91 



-0.20 
-0.32 



257 
340 



M5 



3.82 



14.76 
8.21 



3.41 
2.79 



13.65 
11.26 



2.92 0.034 
2.61 0.065 



-0.68 
-0.37 



0.13 
0.81 



-0.72 
0.47 



117 
171 



M6 



3.82 



24.96 
22.38 



4.04 
3.89 



16.30 
15.72 



1.67 0.051 
4.10 0.032 



0.80 
-0.48 



0.21 

0.88 



0.56 
0.02 



231 
112 



M7 



2.57 



5.15 
21.06 



3.00 
5.16 



8.59 
13.25 



3.09 
2.30 



0.045 
0.077 



-0.54 
0.97 



-0.79 
0.03 



-0.31 
-0.16 



128 
378 



MS 



3.02 



41.84 
34.18 



5.77 
5.39 



17.65 
16.52 



6.25 
5.88 



0.026 
0.031 



0.75 
0.99 



-0.62 
-0.07 



-0.27 
0.09 



116 
130 



TABLE 2 

Initial Positions and Velocities of Merging Dwarfs 



x y z V x V y V z 

Merger (kpc) (kpc) (kpc) (km/s) (km/s) (km/s) 



0.70 
0.76 

1.73 
1.76 

1.85 
1.46 

0.68 
2.00 

1.74 
3.13 

3.61 
4.02 

5.27 
1.29 

1.21 
1.48 



Most classic dSphs have projected ellipticitics, e = 1 — 
b/a (where b and a denote the projected minor and major 
axis of the stellar distribution, respectively) in the range 



0-1 < e < 0.5 (e.g.. lMatedfl998llMcGaugh fc Wolfll20inl) . 

Therefore, we classify as dSphs only those remnants that 
satisfy e < 0.5. For each line-of-sight, we measured e 
at a distance of 2r 1 / 2 , where r 1 / 2 is the half-light ra- 
dius. Such distances are co nsistent with those adopted i n 
observational studies (e.g.. iMcConnachie fc Irwinl l2006') . 
r i/2 were determined by calculating the radius contain- 
ing half the total number of stars in the surface density 
distribution of each remnant. The range of r 1 / 2 in Ta- 
ble [3] corresponds to the minimum and maximum values 
determined for the lines of sight along the three princi- 
pal axes and reflects the fact that the remnants are not 
spherical. Figure [1] shows the surface density maps of all 
eight merger remnants. 

Columns 3-5 of Table |3] list the range of values of 
V^-ot/c* and e for each remnant together with the out- 
come of the classification scheme based on the two cri- 
teria above. Any random line-of-sight would correspond 
to values of V ro t/<J* and e within the ranges quoted in 
Tabled 

Out of eight initial merger remnants, three (M3, M6, 
and M8) would be classified as bona fide dSphs. If we 
slightly relax the shape criterion above by considering 
the most elon gated classic dSph Ursa Minor with e « 0.6 
(|Mateolll99ct L then seven of our remnants would qualify 
as dSphs (we indicate this by using the notation " dSph?" 
in column 5 of Table [3] to characterize remnants Ml, 
M2, M5, and M7). As in mer gers between massive disk 
galaxies (e.g. JCox et al.ll2006h . the values of V m t/o* are 
affected by both the mass ratios of the merging systems 
and the degree of initial alignment of the disks' angular 
momenta. Indeed, despite the near alignment of the disk 
spins in M8, a dSph does form. This is because the disks 



Ml 



-16.69 
18.15 



-17.68 
19.22 



-4.84 
5.26 



4.84 
-5.26 



1.99 
-2.16 



M2 



-6.67 
6.78 



0.29 
-0.30 



-3.13 
3.19 



5.52 
-5.63 



2.36 
-2.40 



M3 



4.78 
-3.76 



-0.04 
0.03 



-0.23 
0.18 



-0.64 
0.50 



-4.18 
3.29 



M4 



-1.96 
5.74 



-3.30 
9.65 



0.54 
-1.58 



-3.04 



0.19 
-0.56 



M5 



-0.02 
0.04 



2.06 
-3.70 



0.18 
-0.32 



-1.10 
1.97 



-6.38 
11.47 



M6 



4.58 
-5.11 



6.34 
-7.07 



3.75 
-4.18 



-3.33 
3.71 



1.99 
-2.21 



M7 



1.49 
-0.36 



2.23 
-0.57 



6.13 
-1.50 



-0.54 
0.13 



-20.50 
5.01 



MS 



-7.76 
9.50 



-3.02 
3.70 



-5.87 
7.18 



4.54 
-5.56 



-6.91 
8.46 
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Fig. 1. — Surface number density distribution of stars in the merger remnants. Results are presented for the most non-spherical views 
along the intermediate axis of the stellar distribution, y (where x and z denote the major and minor axes, respectively). The stars were 
binned into 0.2 kpc X 0.2 kpc fields perpendicular to the line-of-sight. The contours correspond to the number of stars N within such a bin 
and are equally spaced by 0.2 in log N. The innermost contours are in the range logN = 3.2 — 4.2 depending on the remnant. 



TABLE 3 
Properties of Merger Remnants 



T-y/l My [ty Color 

Merger (kpc) Vrot/a* e = l — b/a Classification (mag) (mag arcsec -2 ) in Figure [2] 

(1) (2) (3) (4) (5) (6) (7) (8) 



Ml 


0.22-0.40 


0.10-0.20 


0.11-0.63 


dSph? 


-(8.3-9.2) 


25.1-27.1 


green 


M2 


0.22-0.38 


0.09-0.69 


0.04-0.59 


dSph? 


-(8.5-9.4) 


24.9-26.6 


red 


M3 


0.24-0.32 


0.12-0.70 


0.14-0.45 


dSph 


-(8.4-9.3) 


25.0-26.7 


blue 


Ml 


1.24-1.50 


0.19-1.88 


0.12-0.50 


non-dSph 


-(10.3-11.2) 


26.4-27.6 


black 


M5 


0.30-0.52 


0.20-0.60 


0.03-0.59 


dSph? 


-(9.0-9.9) 


25.1-26.8 


orange 


M6 


0.42-0.60 


0.22-0.75 


0.06-0.50 


dSph 


-(9.7-10.7) 


24.7-26.1 


magenta 


M7 


0.28-0.38 


0.33-0.92 


0.08-0.54 


dSph? 


-(9.1-10.0) 


24.4-25.9 


cyan 


M8 


0.29-0.41 


0.09-0.50 


0.07-0.47 


dSph 


-(10.3-11.2) 


23.7-25.3 


brown 



are effectively destroyed in this nearly equal-mass en- 
counter. On the other hand, mergers M4 and M7 which 
correspond to the smallest mass ratios produce remnants 
that exhibit the highest values of Vrot/c*. In these cases, 
the remnants maintain part of the original rotation of the 
more massive disks which were not completely destroyed 
by the encounters. Moreover, although M7 corresponds 
to a smaller mass ratio compared to M4, its remnant 
exhibits a lower value of V TO t/(J* and is classified as a 
dSph. This is because the intrinsic spins of the disks in 
this case are oriented almost in the opposite direction, 
whereas the disk spins are nearly aligned in M4. 

Columns 6 and 7 of Table[3]list the absolute magnitude 
in the V-band, My, and the central surface brightness, 
fiv- For simplicity, the magnitudes were determined by 
assuming that our imaginary observer will be able to 
detect all stars. The total stellar masses were trans- 
lated into luminosities assuming a stellar mass-to-light 
ratio of M*/Ly = (2.5 ± 1)M q /Lq. The mean value is 



the present-time prediction for a single, low-metallicity 
stellar population i n the standard model described by 
iBruzual fc Cha riot (2003|). The lower limit roughly cor- 
responds to the mean M*/Ly for pressure-supported LG 
dwarfs e stimated fro m resolved stellar population stud- 
ies (e.g.. lMateolll998D . The 1Mq/Lq variation addresses 
some of the uncertainty regarding our assumptions for 
the mean M*/Ly above and the lack of relevant physical 
processes in the simulations (e.g., star formation) and il- 
lustrates how the predicted magnitudes may vary due to 
unknown M^/Ly. 

Lastly, fiy were computed by considering stars within 
0.2 7^/2 from the center of each remnant. Such distances 
roughly correspond to the innermost data points of sur- 
face brightness distribution s of observed dwarfs (e.g., 
IMcConnachie fc Irwin! [2006h . The ranges of /iy in col- 
umn 7 of Table [3] include both the variation of ri/ 2 and 
the uncertainty in M*/Ly. 
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Fig. 2. — Absolute magnitude, My, versus central surface bright- 
ness, fiy (upper panel) and My versus projected half-light radius, 
r i/2 (bottom panel) for dwarf galaxies. The merger remnants 
correspond to color rectangles whose sizes indicate the ranges of 
values in Table [3] Symbols show results for observed LG dwarfs: 
open symbols correspond to dlrrs, filled black symbols show results 
for both dSph and dwarf spheroidal/dwarf elliptical (dSph/dE) 
systems, and filled gray symbols correspond to "transition-type" 
(dlrr/dSph) dwarfs. Isolated dSphs Cetus and Tucana are addi- 
tionally marked wit h red circles. All da ta for observed dwarfs arc 
taken from Table 2 of[Lokas ct al. (2011). The remnants that would 
be classified as dSphs and the LG classic dSphs occupy essentially 
the same regions on both planes. 

4. DISCUSSION 

Galaxies can be characterized by correlations between 
their fundamental parameters. This also applies to LG 
dwarf galaxies which occupy charac teristic regions in th e 
My — fiy an d My — ri/2 planes (e.g.. lTolstov et ail 2009). 
Figure [2] illustrates how our merger remnants (colored 
rectangles) compare with the population of LG dwarf 
galaxies (symbols) in these two planes. Except M4 which 
is not classified as a dSph, our merger products and the 
LG classic dSphs (filled black symbols) occupy essentially 
the same regions on the My — fiy and My — r 1 / 2 planes. 
This further suggests that merging between rotationally- 



supported dwarfs constitutes a viable mechanism for the 
formation of dSphs. 

Using the subhalo catalogues we have traced the his- 
tory of the merger remnants in the cosmological simu- 
lation to z = 0. From the seven systems that could be 
classified as dSphs, four (M2, M3, M6, and M7) sur- 
vived as individual entities to the present-time. Two of 
the remaining merger products (Ml and M8) were ac- 
creted by a factor of ~ 100 more massive halos, while 
M5 participated in a 1:5 interaction with a more massive 
object. Interestingly, all four merger remnants that have 
survived as individual entities to z = appear to be re- 
mote, located in the periphery of the simulated LG at 
distances > 800 kpc from either the MW or M31. These 
systems constitute plausible counterparts of the oddly 
isolated dSphs Cetus and Tucana which reside in the LG 
outskirts, t he former at a distance of ~ 800 kpc from the 
MW (e.g.. iMcConnachie etld1l2005f) while the la tter at 
- 1300 kpc from M31 (e.g.. iSaviane et al.lH996l ). The 
observational parameters of Tucana (My = —9.5 mag, 
ri/2 = 274 pc, fiy = 25.1 mag arcsec -2 , Kot/cr* ~ 1, and 
e = 0.48; Table 2 of lLokas et all 127)11 are in reasonable 
agreement with those of dSphs M3 and M8 in Table [3f. 
On the other hand, the Cetus dwarf (My = —11.3 mag, 
r i/2 = 600 pc, fiy — 25.0 mag arcsec -2 , Kot/c* = 0.45, 
and e = 0.33; Table 2 of lLokas et all 120 111) is brighter 
and more extended and therefore more akin to our dSph 
M6. Overall, the agreement is particularly noteworthy 
because our simulation program did not aim to repro- 
duce the properties of these objects. 

The evolutionary scenario proposed here may also 
be relevant for the relatively isolated trans itional 
(dlrr /dSph) dwarfs Phoenix, LGS3, and Pegasus (jMateol 
1998) that usually differ from dSphs only in the presence 
of small amounts of g as, and even for the re mote peculiar 
dwarf galaxy VV124 (|Kopylov et al.ll2008D . 

Alternative models exist for explaining the puzzling 
presence of Tucana and Cetus in the outskirts of the 
LG. For example, as a result of three-body interac- 
tions, satellites can acquire extremely energetic orbits 
with apocenters beyond the R v \ r of the primary and be 
ejecte d to large distances ()Sales et alJl2007tlLudlow et al.l 
2009). In this model, ejected subhalos are typically the 
least massive members of a group of satellites that is 
tidally disrupted by the host. Tidal interactions with the 
more massive companions within the groups may already 
transform the eventually ejected dwarfs into dSphs, via 
either tidal stirring, p erhaps enhanced by resonances 
(D'Onghi a et al.l [20 09), or mergers. Ongoing investi- 
gations o f the stellar components of Tucana and Ce- 
tus (e.g., iBernard et all 120091: iMonelli et aHl2010al). in- 
cludin g their star formation histories (e.g.. lMoneili etaLI 
l2010bf ) , may soon provide constraints on the various com- 
peting models for their origin. 

Lastly, our simulations neglect gas dissipation and 
star formation. Dwarf galaxies are known to be ga s 
rich at least at low redshift (e.g., iGeha et ail [2006). 
Although such dissipative processes should be impor- 
tant in en counters of massive ga laxies with cold gas 
disks (e.g.. iKazantzidis et al.l[2005l ). they should have a 

7 Due to the very large distances involved the detection of rota- 
tion at such high levels for Tucana should be regarded as tentative 
and needs to be confirmed by future observations. 
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rather marginal effect in our case. Indeed, the dwarf 
halos considered here (Table [1} have both low V v i T < 
20 kms -1 and concentrations (1 < c < 6.5), imply- 
ing maximum circular veloci ties Knax only slightly above 
V v ir- iKaufman n et al.l (|2007l) have shown that halos with 
V max < 35 kms -1 would not host thin gaseous disks but 
rather spheroidal gas distributions. 

Furthermore, neglecting self- shielding, Mayer et al. 
(2006) have shown that the gas in dwarf halos would 
still be heated to T > 10 4 K by the cosmic ionizing 
background at z > 2, when all of our mergers occur. 
Given that K,ir ~ 20 krns" 1 corresponds to T v ; r ~ 10 4 K, 
the gas distribution would form an extended pressure- 
supported envelope. Depending on halo mass, this en- 
velope would be barely confined by the halo poten- 
tial, or even be comp letely evaporated from the halo 
(|Barkana fc Loeblll999D . Regardless, the gas component 
would not remain confined into the central disk, except 
perhaps near the very center where self-shielding might 
be sufficient to preserve a tiny fraction of cold dense gas 
(jSusa fc Umemural 12004) . Therefore, either the merg- 
ers would occur essentially in a dissipationless regime, 
as we have assumed here, or they would contain a tenu- 
ous hot atmosphere which could be shock heated further 
as a result of the merger (e.g., [Kazantz idis et al.l I2005D 
and become unbound thereafter in the absence of effi- 
cient cooling. In either case, no gas inflows and no ap- 
preciable triggered star formation would occur, contrary 



to what happens in merg ers of massive galaxies (e.g,. 
iBarnes &; Hernquistl 119961 ). In summary, the effects of 
gas and star formation on the structure of the merger 
remnants are expected to be weak. Therefore, although 
they will have to be confirmed with hydrodynamical sim- 
ulations, the qualitative arguments outlined above sug- 
gest that our collisionless experiments have likely cap- 
tured the essence of the morphological evolution in merg- 
ers of disky dwarf galaxies. 
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